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Abstract: We examine in detail the relationship between smooth fast quantum 
quenches, characterized by a time scale and instantaneous quenches, within the 
framework of exactly solvable mass quenches in free scalar held theory. Our earlier 
studies [1, 2] highlighted that the two protocols remain distinct in the limit —)■ 0 
because of the relation of the quench rate to the UV cut-off, i.e., l/5t -C A always 
holds in the fast smooth quenches while l/5t ~ A for instantaneous quenches. Here 
we study UV hnite quantities like correlators at hnite spatial distances and the ex¬ 
cess energy produced above the hnal ground state energy. We show that at late times 
and large distances (compared to the quench time scale) the smooth quench correlator 
approaches that for the instantaneous quench. At early times, we hnd that for small 
spatial separation and small ht, the correlator scales universally with ht, exactly as 
in the scaling of renormalized one point functions found in earlier work. At larger 
separation, the dependence on 5t drops out. The excess energy density is hnite (for 
hnite m5f) and scales in a universal fashion for all d. However, the scaling behaviour 
produces a divergent result in the limit m5t —)■ 0 for d > 4, just as in an instantaneous 
quench, where it is UV divergent for d > 4. We argue that similar results hold for ar¬ 
bitrary interacting theories: the excess energy density produced is expected to diverge 
for scaling dimensions A > d/2. 
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1 Introduction 


Universal scaling behavior in systems undergoing a quantum (or thermal) quench which 
involves critical points have been a subject of great interest in recent years [3-6]. The 
classic example of such behavior is Kibble-Zurek scaling [4, 5] which involves a quench 
which starts from a gapped phase at a rate which is slow compared to the scale set by 
the initial gap. At the other extreme, there are a different set of universal behaviors 
in two-dimensional held theories which are quenched instantaneously from a gapped 
phase to a critical point [7, 8] and for instantaneous quenches which can be treated 
perturbatively [9]. 

The AdS/CFT correspondence has yielded important insight in this area, both 
for Kibble-Zurek scaling [10] and for novel non-equilibrium phases [11]. Perhaps more 
signihcantly, holographic studies have led to the discovery of new scaling behavior for 
smooth quenches which are fast compared to the physical mass scales, but slow compared 
to the UV scale [12, 13]. In [1] and [2] we argued that this scaling law holds regardless of 
holography, and is valid for time dependent relevant deformations of generic conformal 
held theories (see also [14]). Consider an action 

S = Soft + j d^x X{t) (Pa(^U) (1-1) 

where the conformal dimension of the operator (P is A and X{t) interpolates between 
the constant values Ai and A 2 (with an amplitude variation of 6X) over a time of order 
6t. Then in the fast quench limit 

St < (1.2) 


the renormalized energy density SSren scales as 

SX‘^ 

Similarly, the peak of the renormalized expectation value of the quenched operator, 
measured at times earlier than or soon after the end of the quench, was also found to 
scale as 

{OA)ren ~ §f2A-d ’ 

This general result emerged out of detailed investigations of exactly solvable mass 
quenches in free bosonic and fermionic theories. One important outcome of our analysis 
was an understanding of the relationship between smooth fast quenches for small St 
and the instantaneous quenches of e.g., [7, 8]. The latter involve a quench rate which 
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is fast compared to all scales, including the UV cutoff, while smooth quench rates are 
faster than any physical scales, but slower than the cutoff scale. On the other hand, 
local quantities like the energy density or {O) involve a sum over all momenta all the 
way to the cutoff — for such quantities one does not expect the smooth quench result 
to agree with those in the instantaneous quench. By the same token, one would expect 
that for correlators at hnite separations larger than 6t, there should be agreement. In 
[1, 2] we also explored if the late time behavior of local quantities also agree, finding 
agreement at least in the d = 3 case. 

In this paper we continue to explore the relationship between fast but smooth 
quenches and instantaneous quenches in further detail. Our analysis will focus on 
quenches in free scalar field theory with a time-dependent mass. However we argue 
that the lessons we draw there will be valid for quenches in interacting theories of the 
type described above. 

The scaling laws in [1, 2] were derived for renormalized composite operators, which 
are the appropriate quantities for quench rates much slower than the cutoff scale. In 
this work, we examine the late time behavior of such operators. In addition, we focus 
on quantities which are UV hnite, e.g., two-point correlation functions at hnite spatial 
separations and the excess energy over the ground state energy at late times. 

In section 3, we consider late time correlators, t 3> St. We will show that for (suit¬ 
ably dehned) large spatial separations, these correlators agree with the correlators for 
an instantaneous quench. For separations r which are very small, i.e., mr -C 1 there is 
once again agreement, rehecting the fact that the dominant singular behavior for small 
separation is independent of any time dependence of the mass. The corrections to this 
leading small distance behavior are in one-to-one correspondence with the countert¬ 
erms necessary to renormalize the composite operator 0^. In particular, the sub leading 
small distance divergences involve derivatives of the mass function for d > 6. For 
intermediate separations, the two quench protocols lead to genuinely different results. 

In section 4, we turn our attention to correlators at finite times t ^ St and show that 
for rSt ^ 1 the correlator becomes independent of St as expected. For m~^ > St > r we 
hnd that the correlator exhibits a scaling behavior identical to that of the renormalized 
local operator (0^). 

In section 5, we consider the renormalized local quantity ((/)^) at late times. We 
show that this quantity agrees for both quench protocols only for d = 3. For d = 5 and 
hnite St, there is a slight difference between the smooth and the instantaneous answer 
in the limit of St —)■ 0, while for the instantaneous quench, (0^) is UV divergent for 
d > 5. 

In section 6, we consider the difference of the energy at late times and the ground 
state energy with the hnal value of the mass. This is one measure of the excess energy 
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produced during the quench. We show that this quantity is explicitly UV hnite. For 


d < 3 this becomes independent of 6t, in the m6t -C 1 limit. The next order correction, 
which scales as a power of 6t, is identical to the behavior of the renormalized energy 
in [1, 2]. For d > 4, the energy diverges in the dt —)■ 0 limit, in the same way as the 
renormalized energy considered in [1, 2]. 

In section 7, we discuss the validity of our results for the excess energy produced 
for general interacting held theories. 

In section 8, we conclude with a brief discussion of our results and also consider 
various possible measures of the energy produced by the quench and their relationship. 

2 Bogoliubov coefficients for smooth and instantaneous quenches 

Consider a scalar held in d space-time dimensions with a time dependent mass. 



( 2 . 1 ) 


This theory is exactly solvable for a variety of diherent mass prohles, as described in 


[1, 2]. The quench protocol which we focus on here, involves the mass going from an 
initial value m to zero at late times over a time scale 6t with the smooth prohle 



( 2 . 2 ) 


To solve the Klein-Gordon equation, we decompose the scalar held into momentum 
modes 




The exact solution of the held equation is given by [1, 2, 15] 



(2.4) 


(2.5) 


where 2 F 1 is fhe usual hypergeometric function and 



(2.6) 
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The modes are the “in” modes: they behave as plane waves in the inhnite past and 
the ag annihilate the in-vacuum, i.e., ag|m,0) = 0. There is also another set of modes, 
the “out-modes”, which become plane waves in the inhnite future, 

= , exp(ifc • X — iu:+t — iu^St log(2 cosht/^t)) x 

V‘^^out 

2 F 1 -h iu-6t, iuj-6t] 1 -f iUoutSt] -—^ _ ^ 2 . 7 ) 

In terms of these, the held operator is 

^ = J (2^)(d-b/2 {h > ^here [bj:, &[-] = 5'^-\k - k '). (2.8) 

The Bogoliubov transformation that relates these two sets of modes is given by [15] 


where 


= 






k k k ^k —k'> 


^OUt r(^l 

^in T{—iu+6t)T{l — iuj+6t) ’ 
^out r(f 

UJi: 


The Heisenberg-picture state of the system is the “in” vacuum. 


(2.9) 


( 2 . 10 ) 


a^|m, 0) = 0 . 


( 2 . 11 ) 


We will be interested in analysing several quantities: (i) the two-point correlator of the 
held at a hnite spatial separation (ii) the expectation value of the composite operator 
(jF and (iii) the net energy density produced. In fact, the rate of change of the energy 
density is related to (0^) by the Ward identity 

dt{8) = ]^dtrri^{t){(t)^). (2.12) 

The two point correlation function of the scalar held under the quench reads 


(m, O|0(T, t)0(x', T)|m, 0) = 


r d8 ^k ^ 

J ( 2)^1 


r d'^-^k 

J (27r)'^-i 




(2.13) 
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We will be interested in the relationship of the results of a smooth quench as in 
(2.2) to that of an instantaneous quench from a mass m to zero mass, 

'^instant(t') Tn9(^ t). (2.14) 


The “in” and “out” modes for such an instantaneous quench have a trivial plane wave 
prohle for before and after the quench, respectively 

instant 

instant 


V2 




'sJ^U^out 


^ t < 0 , 

I 

^i(k-x-uJoutt) ^ > g 


(2.15) 


The Bogoliubov coefficients for the instantaneous quench are determined by demanding 
that the mode functions and their first derivatives are continuous at t = 0. This yields: 


a 


instant 


UJ+ 


OJin^out 


and 




instant 

k 


OJ- 


y/OJin^out 


(2.16) 


The correlator for an instantaneous quench can be easily computed using (2.16) (or by 
directly matching the operator solutions across f = 0 as in [8]), 


^ ^ -IXinf 

‘ZUout 


^ik-{x—x') 

COS U out {t - t') + 


(2.17) 
eosa;,,,(t + T) . 

J 


The comparison of the instantaneous quench with the smooth fast quench is only mean¬ 
ingful at late times when the variation of the mass in the smooth quench is over, i.e., 
when t S> and t' S> 5t. For such times, the mode functions vj:{x, t) —)■ which 

are exactly the mode functions for f > 0 in instantaneous quench. 

In what follows, it is useful to look at the behavior of the Bogoliubov coefficients 
in various regimes. 


1. First consider the limit 

1 • 

It may be easily checked that in this limit, the smooth quench Bogoliubov coef¬ 
ficients (2.10) reduce to the instantaneous quench coefficients (2.16), regardless 
of the value of \k\/m. This means that the smooth quench approaches an in¬ 
stantaneous quench when 6t is small compared to all other length scales in the 
problem. In particular, this means that the momentum space correlators at some 
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momentum k will approach the instantaneous correlator only when m6t 1 as 
well as -C 1. In [1] we discussed the implications of this for expectation val¬ 
ues of local quantities like or the energy density. Generically, the small 

m6t 1 limit of these local quantities will not agree with the instantaneous 
quench result since these quantities involve an integration over momenta all way 
upto the cutoff, and the physical smooth quenches, in which we are interested, are 
fast compared to physical mass scales but slow compared to the cutoff scale. It 
can be the case, however, that the integrand is a rapidly decaying function. If so, 
even if we have to integrate to arbitrarily high momenta, the main contributions 
will come from low momenta and then there would be agreement between the 
two protocols. For example, as shown in [2], this is what happens in evaluating 
(0^) in d = 3. In section 5, we will go back to this discussion and show that in 
higher dimensions this is not generally true. 


2. Now consider the limit 

mdt <C 1, \k\/m^l, I= arbitrary. (2.19) 


Once again in this limit, the Bogololiubov coefficients (2.10) approach the instan¬ 
taneous quench coefficients (2.16), for any hnite value of \k\6 t. In fact, in this 
limit the coefficients (2.10) behave as 

9 

Tfl 

—)■ 1 — [1 ~ ik6t'rp{l — ik6t) + ik6t'ip{—ik6t)] + 0{m'^/k ‘^), 

9 

771 

r(*fc^t) + 0{m^/k ^), (2.20) 

where 'ip{x) denotes the digamma function, i.e., ijj{x) = cIa;logr(x). For the 
instantaneous quench, instead, they behave as 


instant 


->■ 1 - 1 - 0{m^/k ‘^), 


/S 


instant 


-)■ 


2 

m 

4P ■ 


( 2 . 21 ) 


Thus to leading order in m‘^/k‘^, we have ag = l,/3g = 0 for both smooth and 
instantaneous quenches regardless of the value of k6t. This is a reflection of the 
fact that very high momentum modes are not excited by the quench, i.e., to 
leading order the quench is immaterial for these modes. The subleading terms in 
(2.20) are of course dependent on 6 t. In fact, for hnite k6t, the subleading term 
in is 0{m?/k"^). However, if in addition we have k5t 1, this 0{m?/k"^) term 
is cancelled, as it should be. 


In the next section we discuss the implications of these observations for real space 
correlation functions. 
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3 Late time spatial correlators 

In this section, we examine equal time correlation functions at finite spatial separations 


C{t,f) = (0(f,f)0(f,O)), 


(3.1) 


and compare the result for smooth fast quenches and instantaneous quenches at late 
times. For simplicity, we only explicitly consider the correlators in odd spacetime 
dimensions in the following. We will consider this correlator in eq. (3.1) in three 
different situations: the hrst one is the equal time correlator for a smooth quench from 
a mass m to zero mass, as in eq. (2.2), 




smooth 


2(27r 


\d—l 


d^-^k 


Ak-r 


ar 


+ lA 


.1 


+ e 


—2ikt 1 


(3,2) 


where and fdj: are given by eq. (2.10). The second quantity is the equal time correlator 
for an instantaneous quench which can be read off from eq. (2.17), 

C(t, r).„„ = J <f-^k e-‘-' (p + sm2(it))^ , (3,3) 

This correlator was studied in e.g., [8]. Finally, we consider the correlator for a constant 
mass m = 0, 

Cconsti^ = , , f (3.4) 

constK J 2(27r)'^-i J \k\ 

Constant mass correlators are evaluated in detail in Appendix A, including the case of 
m = 0 — see eq. (A.8). Performing the angular integrals above we find 


C{tt smooth / t^smooth(^k,t,T^ dk 


Or. 


1 


d-3 

arT 2 


dk k 2 J d-3 (kr ) -j |q(^| + |/ 3 ^| + e * + e 


(3.5) 


2 I ^ /QA 2 ikt j_ —2ikt 1 


C{tiT^instant / ^instanti^k^t^T^ dk 

J 

= -^ j dk k~ Jd-3 (kr) 


(3.6) 


k^ + m? sin^(/cf) 


Or.r 2 


Cconsti^^') / tlconstik')'^'^ dk 


2 ' ' k\/k‘^ + m? 

1 


d-3 


Or 


d-3 

OrT 2 


dk k ‘2 J d-3 {kr ), 


(3.7) 


where cTc = 2 2 tt 2 and J d -3 is the Bessel function of order 

2 ^ 













Before proceeding with the detailed calculations, let us present our intuitive ex¬ 
pectations for the comparison, as well as the results found below. As in section 2, we 
are considering quenches which take the mass from some hxed initial value m to zero 
after the quench. There are several different dimensionful parameters at play in our 
correlators, i.e., t, 6t, r = Ir*! and m, and so hrst, we wish to clearly specify how the 
corresponding scales are ordered in our considerations below. First, for the smooth 
quenches, we are considering the fast quench limit and hence we have m6t 1. We 
are also examining the correlators in late time limit and hence 6 t ^ t. While these 
inequalities do not £x a relation between m and t, we will only present results for the 
case mt > 1. That is, the following discussion explicitly considers quenches where 

6 t <C 1/m < t. (3.8) 

We have verihed that the general behaviour is the same in regimes where mt < 1, as 
long as the inequalities for the fast quench and late time limits are both satished. 

Given the ordering in eq. (3.8), we still have one remaining scale unspecihed, namely 
the spatial separation r = |r |. In the following, we compare the correlator (3.1) for fast 
smooth quenches and that (3.3) for instantaneous quenches for r in different regimes. 
The natural intuition is that in the Fourier transform, only momenta satisfying fc < 1/r 
will contribute signihcantly to the correlator. Hence given that we are in the fast quench 
regime with m6t 1, if we further choose r ^ then both of the inequalities in 
eq. (2.18) should be satished for the momenta contributing to the correlators. Hence 
the analysis in the previous section would indicate that at late times, the integrand in 
eq. (2.13) for the smooth quenches effectively reduces to that in eq. (2.17) matching the 
correlator for an instantaneous quench. We will explicitly verify that this expectation 
is realized by numerically comparing the full expression (2.13) for smooth quenches 
with the instantaneous quench result (2.17) at late times. 

Continuing with the above intuition, differences between the late-time correlators 
for the two different quench protocols are expected to arise in a regime where the spatial 
separation is comparable to the quench time, i.e., r < St. In view of eq. (3.8), this means 
that we would be examining the correlator at short distance scales. However, we also 
found that for very large fc 3> m the leading behavior of the Bogoliubov coefficients 
are in fact independent of the quench rate — see discussion following eq. (2.19). 
This immediately implies that the leading singularity at small r is independent of any 
quench protocol and therefore one gets the same singular behaviour as the constant 
mass correlators in this regime, namely, C{t,f) cx — see Appendix A. Hence 

the leading behaviour in the correlators produced by the smooth and instantaneous 
quenches again agrees in this regime. As we will show below, the subleading singularities 
are in fact different in high dimensions. 
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The interesting regime where the difference between a smooth and instantaneous 
quench would show up is therefore the intermediate values of r. To study this differ¬ 
ence it is convenient to eliminate the leading small-r contribution by calculating the 
difference between the two quench correlators, i.e., Csmoothit,f) — Cinstant or al¬ 
ternatively by subtracting the fixed mass correlator from each of the quench correlators 
individually. As we describe in detail below, this difference of the late-time correlators 
indicates that the subleading behaviour, in fact, also agrees for d = 3 but a small finite 
difference arises for d = 5. For d = 7 and higher dimensions, the difference still diverges 
in the limit r —)• 0. 

3.1 Numerical results for various dimensions 

We now evaluate the k integrals in eqs. (3.5) - (3.7) numerically. However the integrands 
are typically oscillating rapidly with a growing envelope. Hence to evaluate the integral 
properly, we need to regulate the integral. We do so by introducing a convergence factor 
exp(—a/c) with a > 0 and evaluating the integral in the limit a —)■ 0, as discussed in 
Appendix A. Using this regulator, the integral in eq. (3.7) can be evaluated to yield 
Cconsti'k) oc as in eq. (A. 8). This is the standard massless propagator in d 

dimensions. For our numerical calculations, we typically use a = 10“^, which we can 
verify is small enough that the integrals accurately converge to their limiting values. 

As noted in the above discussion, for large k ^ m, Csmoothik) and Constant {k) 
are essentially identical to Cconstik), implying that the leading divergence in all the 
corresponding correlators is for small r. Hence, the integrands are very close to 

each other for a large range of /c, as illustrated in fig. 1. 

Therefore, in order to highlight the differences between the smooth and instanta¬ 
neous quenches, we will subtract off eq. (3.7) to define for both cases 

C{t,r) = — [ C{k,t,r)dk = C{t,r) - Cconst{r) = — [ {C{k,t,r) - Cconst{k,t,r)) dk . 

(3.9) 

In terms of the integrands, subtracting Cconstik) suppresses the growing oscillations at 
large k that, e.g., we see in fig. 1. In the full correlator, this subtraction removes the 
divergent behaviour at r —)■ 0, which makes it easier to identify differences in the finite 
remainder. If this behaviour was not removed, for instance, both the instantaneous and 
the smooth quench would both grow as in exactly the same way as r —>■ 0 and 

it would be extremely difficult to find any differences between the correlators for the 
two different quench protocols in this regime.^ In comparing the subtracted integrands 

^This would be analogous to analysing the bare expectation value instead of the renormalized 
quantity (d^)ren in [li 2]. Of course, the interesting physical quantity is in the renormalized expectation 
value — see further discussion in section 3.2. 
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Figure 1. (Colour online) Integrands Cinstant (blue), Csmooth (yellow) and Cconst (green) of 
the instantaneous quench, smooth quench and constant massless case, for d = 5, t = 10, 
r = 10 and = 1/20 as a function of momentum k — all dimensionful quantities are given 
in units of the initial mass m. All three integrands are nearly identical for large momenta. 
Further the instantaneous and smooth quench curves overlap at all scales. 

C below, we start by considering d = 5. 

In fig. 1, as well as the good agreement between the qnench correlators and the 
constant massless correlator at large k, we see that the two integrands, Cgmoothik) 
and Cinstant{k), agree for thronghout the momentum range shown there. Given our 
discussion at the opening of this section, we only expect differences to arise when 
k >1/61. One way to illustrate these differences is to make St larger, as illustrated in 
£g. 2. Panel (a) shows the subtracted integrands in the range 0 < A; < 2 with the same 
parameters as used in fig. 1, i.e., mt = 10, mr = 10 and mSt = 1/20 for d = 5, and the 
two curves precisely overlap in this momentum range. The only change of parameters 
in panel (b) is that here mSt = 1/2 and we clearly see that small differences appear 
between the integrand for the smooth quench and that for the instantaneous quench. 
Note, however, that with m6t = 1/2, the smooth quench only barely satishes eq. (2.18) 
and so while useful to illustrate possible differences, this example is not really in the 
fast quench regime. 
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(a) St = 1/20 


(b) St = 1/2 


Figure 2. (Colour online) Subtracted integrands as a function of (small) momentum k. In 
this case we are plotting for d = 5, t = 10, r = 10 (with the units set by m). The yellow line 
corresponds to the instantaneous quench while the blue one to the smooth. Panel (a) shows 
that no detectable differences appear with 6t = 1/20 in the range 0 < k < 2. However, in 
panel (b), minor differences occur in this range when 6t = 1/2. 


Focusing on the parameters mt = 10, mr = 10 and mSt = 1/20 (for d = 5), we 
see in fig. 3 that the subtracted integrands continue to agree for much larger values 
of k. However, with k/m ~ 20, differences are clearly visible in panel (b). However, 
comparing the vertical scale in these two plots to that in panel (a) of £g. 1, we see that 
these large k contributions to the subtracted integrand are highly suppressed.^ Hence 
the differences should not be expected to contribute to the full integral, i.e., they should 
not produce signihcant differences in the position-space correlators. 

^Further, these large k contributions are also highly oscillatory, so they will tend to cancel out in 
the integral. 
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(a) Intermediate momentum 


(b) Large momentum 


Figure 3. (Colour online) Subtracted integrands as a function of momentum k. In this 
case we are plotting for d = 5,t = 10,r = 10 and 5t = 1/20 (with the units set by m). 
The blue line corresponds to the smooth quench while the yellow one to the instantaneous 
quench. Panel (a) shows an intermediate regime where no significant differences between the 
two integrands are apparent. Visible differences appear for larger k > l/<5t, in panel (b). 



Figure 4. (Colour online) Regulated integrands as a function of momentum k for d = 5 and 
small separation. The blue line corresponds to the smooth quench while the yellow one to the 
instantaneous. At small distances the two integrands are clearly different from each other, 
even at scales k6t < 1. 


It turns out that examining the subtracted integrands for d = 3 yields essentially 
the same results. On the other hand, the situation is also similar in d = 7 for long 
distances. As we decrease r, for d = 7, we see analogous effects to those in £g. 4, i.e., 
the two integrands become clearly different. However, as we will see, after integration, 
the behaviour of the correlator at small distances is very different depending on the 
dimension. 
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100 



m r 

Figure 5. (Colour online) Difference between the late-time correlators for smooth and in¬ 
stantaneous quenches as a function of the separation distance r. The blue line corresponds 
to the d = 3 case while the yellow one belongs to d = 5 and d = 7 is shown in green. We are 
using mt = 10 with mSt = 1/20. In d = 3 and d = 5, the difference remains small for any 
value of r, while in d = 7, it seems to diverge as r —)• 0. 


Given that we understand the behaviour of the (subtracted) integrands, let us 
now compare the (position-space) correlators for instantaneous and smooth quenches 
at different values of r, as shown in fig. 5. In the figure, we have chosen mt = 10 and 
m6t = 1/20 while mr varies from 10 to 0.001. As we have expected, the hgure shows 
that the difference between the correlators goes to zero, or at least is of order 0(5f^), at 
large mr. The results for d = 3 indicate that the difference remains vanishingly small 
for all values of r. In the case of d = 5, a small but hnite difference appears to develop 
as r goes to zero. The differences in £g. 5 are most evident for d = 7. In this case, 
the relative difference is already of order one at mr = 0.1 and the trend in the hgure is 
that it continues to grow at smaller r. Our expectation is that in fact, this difference 
will in fact diverge as r —)■ 0. This conclusion comes from comparing to the late-time 
behaviour of expectation value {4>‘^)ren in fhe next section. 

In summary then, we have explicitly shown that at times long after the quench, the 
correlators generated by instantaneous and the smooth fast quenches are identical at 
large separations. As might be expected, differences only appear at separations of the 
order the quench time 6 t. Further these differences are small in lower dimensions, e.g., 
d = 3,5, but can be quite dramatic in higher dimensions. Interestingly, for d = 3 and 
5, the subtracted correlator of both the instantaneous and the smooth quench (with 
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the mr = 0 piece subtracted out) agree to a rather high precision for any distance mr. 
As we increase mr both answers become even closer. 


3.2 Small r behavior and counterterms for local operators 


In [1, 2], we studied the UV divergences appearing in (0^) with a momentum cut-off 
in great detail. Intuitively, one can think that the correlator at small r provides a 
point-splitting regulator of the same quantity. Hence, the divergences in the correlator 
at small r should be related to the UV divergences of the local quantity (0^). In this 
subsection, we make this connection precise. 

Let’s start with the constant massless correlator, (3.7). For small r we can expand 
the Bessel function to get. 


J d-3 ikr) 
2 




O {{krY) 


(3.10) 


Inserting this in eq. (3.7), we get an integral of exactly the same as the leading 
divergence in (0^). 

Let us hrst recall the set of counterterms that we found in [1, 2] for (0^). There the 
leading divergences in momentum space were determined by performing an adiabatic 
expansion in time derivatives and then expanding our results for large momentum. We 
then found that this was the same as expanding the integrand of (0^). So, to be more 
specihc, in the “in”-basis, the integrand took the form k'^~‘^ojlY\‘ 2 Fi\^ and for large k 
we found [2] 


.(i-2 




2Fi\ 


~ fc'’* ^ — ^7— m^(f) -I- (3m‘^(f) -|- (3-11) 

^ O 

Ud-9 

-(l0m®(f) -I- dfm'^{t) + lOm^(f) + 5dtm‘^{t) dtm‘^{t)) 

32 


Above, we are giving all the terms needed to regulate the expectation value up to 
d = 9. Apart from the divergent terms in the constant mass expectation value (that 
will appear as zeroth order in the adiabatic expansion), they include terms with time 
derivatives of the mass prohle which produce divergences in (0^) for d > 6. These 
(perhaps surprising) terms were carefully analysed in [2]. Now we can express the bare 
expectation value of the local operator in terms of an energy cut-off A, obtained by 
integrating the momentum integral up to a maximum value kmax — A. This yields 


( 02 ) ^ 


1 / A'^ 2 A'’* ^ m^{t) A'^ ® 3m^(t) + d^m?{t) 

2rf-i7r^r (^) U-2 ~ d-A 2 ^ d-6 8 

j ^4^2| X0m2(t) d^m?{t) + bdtm?{t) dtm?^) 


+ ■■■ 
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Working in terms of the “in” modes, the smooth quench correlator can be expressed 


as 

I r ^(d-l)/2 

C(t 1 f" ^smooth d —3 I J d—3 (kv') |2-fl| ) 

acr~ J ^in ^ 

1 f 1 

= -^ / dk-^Jd^{kr) -laFil^ , (3.13) 

acT 2 J fc 2 2 (^in 

where in the last line we have presented the integrand in a way which allows us to make 
use of the above expansion in eq. (3.13). So we have all the ingredients to take the 
limit from the spatial correlator to the expectation value. In particular, if we take the 
r —)■ 0, then the Bessel function can be replaced for its zeroth order expansion shown in 
eq. (3.10). The powers of kr in eq. (3.10) will cancel those same powers appearing in 
eq. (3.13), while the numerical factors will turn the cxc into a ag-^ Then we can use eq. 
(3.11) to expand for large k and hnd that the correlator behaves exactly in the same 
way as (0^) when r —0. In fact, 

C{t,f) -)■ — [ dk (k^~^ — - — m?{t) + -— (3m^(t) + d^m‘^{t)) (3.14) 

^^0 as J \ 2 8 ^ 

k^-Q 

—(l0m®(t) + dfrn?{t) + lOm^(t) dfrn^{t) + bdtrn?{t) + ■ ■ ■ 



Of course for any hnite (positive) r the correlator will be UV hnite. So in principle 
we should be able to perform the integral over momenta to obtain a series expansion in 
inverse powers of r. We now show that this small-r expansion of the correlator is directly 
related to the large-fc expansion for the expectation value. Let us take expression in eq. 
(3.13) and replace k'^~‘^uj~^\ 2 Fi\^ with the expansion given in eq. (3.11). The integrand 
is composed of a series of terms which are constants (not functions of momentum) 
multiplying the Bessel function and some power of k. For hxed value of r, we can 
integrate that expression, using 



k°'Jd-3 (kr) dk 
2 


2 “ r( 4 (d + 2 a-l)) 

^ r hid - 2a - 1))' 


(3.15) 


Now, as an example, consider the leading divergence, i.e., k^ which gives 


1 

d-Z 

a^r 2 



k 2 J d-3 {kr) dk 


I 2‘yr(|-1) 1 ^ 

(Jc ' 


(3.16) 


^We remind the reader that ag is a numerical constant that depends on the space-time dimension d 
and was used in [1, 2] to normalize the expectation value of Explicitly, as = ^ ^ , 

where ^d -2 is the angular volume of a unit (d — 2)-dimensional sphere. 


16 

















after using eq. (3.15) with a = {d — ?>)/2 and some algebra. This shows that UV 
divergences of the local quantity (0^) appear as inverse powers of r in the hnite spatial 
separation correlator, i.e., r plays the role of a point-splitting regulator, replacing the 
momentum cut-off A in eq. (3.13). In the above example we showed that the leading 
divergence is proportional to We can proceed to do the integral for a = — 2. 

This would correspond to divergences proportional to k'^~^ and will lead to a term in 
the spatial correlator which is inversely proportional to Note that in general, 

eq. (3.15) maps the integral over fc" to the power l/r"+^. Also note that there is 
an important difference between the leading divergence and the rest of the divergent 
terms. All of the subleading divergences are proportional either to the instantaneous 
mass m{t) or to time derivatives of m(f), while the leading divergence is independent 
of m{t). This means that the leading term as r —)■ 0 will be inversely proportional to 
but then there will be subleading terms inversely proportional to 
that will also contain factors of the mass and its derivatives. Explicitly we hnd 

... 1 rn^jt) 1 3m^(t) + a,W(t) 

’ Svr'^A y 2 3((i — 6)r^~^ 8 

1 10m®(t) -|- dfm?{t) + lOm^(t) d^m?{t) + bdtm?{t) dtm?‘{t) 

~15(d-6)(d-8)r'^-8 32 ^'' 

In particular, note that to analyse the structure of the correlator in d = 7, we would 
have to take into account a term that is proportional to the second derivative of the 
mass, that will increase as 1/r when r —)■ 0. This will be important to understand the 
behaviour of the correlator in section 4. Also note that this term is also subleading 
compared to the term proportional to m^/r^, that comes second in the expansion of 
eq. (3.18). 

Finally let us note that, even though the momentum cut-off expression in eq. (3.13) 
and the r expansion expression in eq. (3.18) are similar in form, there is no simple way 
to relate the energy scale A with the point-splitting regulator 1/r. Rather, equating 
these two equations we get. 


r(d-i)3^ / (i-r(d- 2 )r(d-i) <^- 2 ) 

A ~ - I 1 -|- - 2 - —rn?{t)r^ + 

2(d-4)r(d- 1)^ 


(3.18) 


This simply points-out that these divergent terms are unphysical and that these two 
regularization schemes have slightly different counterterms. 
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4 Universal scaling in quenched spatial correlators 


In [1, 2] we found an interesting set of universal scaling relations for the expectation 
value of the quenched operator (0^) and the energy density. For the quenches considered 
in this work, this scaling takes the form {4>‘^)ren ~ We also found analytic 

leading order expressions for this expectation value in the case of 6t —)■ 0. For odd 
spacetime dimensions, we found 

+ (4.1) 

which reproduces the above scaling with the mass prohle 2.2 we are considering. On the 
other hand, we have already shown how the spatial correlator approaches the expecta¬ 
tion value of as the separation distance goes to zero. Then an interesting question 
to ask is whether there are any signs of the universal scaling in the spatial correlator. 

To investigate on this question here, we will concentrate on early times, since this 
is the regime where scaling of the local quantities hold. Now we need a suitable object 
to compute. We remind the reader that the scalings hold for renormalized expectation 
values. Given that the bare expectation values are UV divergent, we had to add suitable 
counterterms to eliminate those divergences. On the contrary the spatial correlator is 
hnite for any hnite separation r. However, as discussed in the previous section, the 
counterterms are in precise correspondence with the small r expansion of the correlator. 
In particular the leading UV divergence of (0^) goes as which reflects the leading 
small r divergence of the correlator behaving as So from this perspective, the 

scaling behaviour is only exhibited in a higher order term, which remains hnite as r —)■ 0. 
However, we may still see the scaling in the correlator by subtracting a suitable hxed 
mass correlator to remove the terms which diverge in the small r limit. It turns out 
that in order to eliminate these divergences (which are proportional to powers of m?) 
the interesting object to compute is the difference of the spatial quenched correlator at 
time t with the correlator at a constant value of the mass equal to the instantaneous 
mass at that time t. The latter hxed mass correlator has been computed in Appendix 
A and one hnds. 


Cfixed(j"^ 


1 

d-3 

O'er 2 


k'^dk 


\/k‘^ + mP'it) 


Jd-3 {kr) 


(4.2) 


In order to numerically integrate the correlator, we will go back to the “in” basis — 
see eq. (3.13)— to obtain 


C{t,r) 
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Ud-l)/2 

dk 
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(4.3) 
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First consider the correlator at t = 0. In fig. 6, we computed C(t = 0, r) —Cfixedir) 
for a wide range of values of r and 5t for d = 5 and d = 7. We see a very interesting 
behaviour. Basically we can divide the correlator in three different regions: (i) r > 6t, 
(ii) r < 6t < m~^ and (iii) r < 6t ^ m~^. 

For 6t < r we see that the correlator is essentially a constant independent of 6t that 
depends on r. In fact it turns out that this constant value is exactly the same value of 
the instantaneous quench correlator evaluated at t = 0. Recall that the instantaneous 
quench correlator at t = 0 is exactly the same as the constant mass correlator with 
This coincidence might lead the reader to think that this behaviour is 
something special for t = 0. However, in what follows, we will show that is the general 
behaviour of the correlator at any finite t/5t, as long as r ^ 1/m. 

Let’s start by fixing the dimensionless time r = t/5t. Now we want to analyse the 
r and 5t dependence of the following object: 


C{t,f)-Cfixed{f) 


1 

d-Z 

a^r 2 


dkk 2 Jd-3 (kr) 



1 

^/k'^ + m? 



(4.4) 

Note that the first term inside the big brackets has m^„ in the denominator while the 
fixed mass part carries an m^ equal to that at the particular time we are considering. 

The second important thing to notice is where is the time-dependence in the 
quenched correlator. The only place where r appears explicitly is in the last argu¬ 
ment of the hypergeometric function. Recall that 


2F1 


2-Fi 


iu-5t; 1 — iooindt] 


1 -|- tanh(t/5t)\ 

2 )■ 


(4.5) 


Then, if we fix r by inserting any finite number (or zero) in that last argument, we 
are left with an object that depends only on 6t and r, and we can take the desired 
limit without any problem. So consider now St ^ r. We get that limit by considering 
the limit of the hypergeometric function with St —?• 0. As the second argument is 
proportional to St, to lowest order all the terms in the infinite sum will vanish but the 
first one, so we just get that when St —)■ 0, 2 F 1 = 1 + 0{St). Note that this argument 
is valid only in the case where we fix r and effectively there is no St dependence in the 
last argument. Then, after taking this first limit, we are left with 


C{t,f) - Cfixedif) ^ 


1 

d-Z 

acT 2 


dkk Jd-3 (kr) 



1 

\/k^ + m? 



(4.6) 

But this is nothing more than Cfixedijn^ = — Cfixedijn^ = m^(r))! So, at any time 

we get that for St <^r our object becomes the difference of two fixed mass correlators. 
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In particular, it becomes independent of St and that explains the horizontal dashed 
lines of £g. 6. 

However we can even go further and consider the limit of r -C 1/m. As it is just the 
limit for hxed mass correlators we can extract it directly from eq. (A.7) in Appendix 
A. As we are subtracting two hxed mass correlators, the leading term in the expansion, 
i.e., the one proportional to will cancel and the leading contribution will come 
from the hrst subleading term: 






1 

CTc 



m^(r)\ 

J 


(4.7) 


In our plots of fig. 6, we have m^^ = 1 and m(r = 0) = 1/2, so the difference of 
correlators in the limit of r -C 1/m should go as which agrees with 

the values that the horizontal lines take in the plots. We can add that the agreement 
gets more exact as r takes smaller and smaller values. 

Note, however, that eq. (4.6) is valid for any value of m, as long as r/St 3> 1, but 
eq. (4.7) also needs rm 1. If we concentrate in a regime where rm ~ 1 (see bottom 
line in each plot of £g. 6), then we should expect eq.(4.6) to hold rather than eq.(4.7). 
In fact, for r = m = 1, the bottom dashed line in kg. 6 corresponds to approximately 
0.10 and 0.20 for d = 5 and d = 7, respectively. This is in perfect agreement with eq. 
(4.6) and clearly differs from the 0.25 that eq. (4.7) is predicting. 

To further support our claim we provide equivalent plots but at different times for 
d = 5. In kg. 7, we examine the results for r = 1/2 and r = 1/4. In these cases, the 
expectation is that the limiting value would be 


C{t,r) - Cfi^edir^ 


m}^ — m?{T) 1 — (1/2 — l/2tanhr) 




CTcT 


(4.8) 


and that’s what arises in this kgure.^ 


The second regime r < St < 1/m leads to the most interesting behaviour. The 
correlator now exhibits exactly the same scaling as the 0^ expectation value. The solid 
lines we see in the plots are exactly the lines that come from evaluating eq. (4.1) for 
d = 5 and d = 7 with the present mass quench prokle — see eq. (2.2). This means 
that in this regime exactly the same universal scaling we’ve been discussing for is 
reproduced in the spatial correlator. 

Finally, at least for d = 5, we see that when St ^ 1/m this behaviour breaks down 
and our calculation goes away from the universal scaling line. It would be interesting 

^Forr = 1/2, then C{t,r)—Cfixed{r) « 0.365529/((Tc r) and for r = 1/4, the value is 0.31123/(crc r). 
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to understand better this slow regime as it could be connected to other set of important 
universal scalings in quantum quenches: namely, the Kibble-Zurek scaling [4, 5] and this 
would give a connection between the fast and the slow regime universality in quantum 
quenches. We hope to report on this in the near future. 



(a) d = 5 
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(h) d = 7 


Figure 6. (Colour online) Spatial correlator under a smooth quench at t = 0 as a function 
of both 5t and the distance separation r. In each case, we are subtracting the fixed mass 
correlator with = 1/2. The dashed lines correspond to computing the instantaneous 
quench correlator at t = 0 for the different separations r, that is the same as computing the 
fixed mass correlator with m = rriin = 1. The purple solid line shows the analytic leading 
order contribution to ((/>^), given by eq. (4.1). 
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(a) T = 1/4 



r= 1 




r = 

r = 


10^ 

1 

10^ 

1 

10“ 


(b) T = 1/2 


Figure 7. (Colour online) Spatial correlator under a smooth quench for hxed t/5t = r as a 
function of both 5t and the distance separation r. In each case, we are subtracting the fixed 
mass correlator with m?{t/5t = r). The dashed lines correspond to computing the fixed mass 
quench correlator with m = rriin for the different separations r. The purple solid line shows 
the analytic leading order contribution to (0^), given by eq. (4.1). 


5 Late time behaviour of 

In [2], we found some interesting behaviour for the expectation value of 0^ when we 
examined d = 4 at late times. Essentially, the expectation value for the smooth quench 
did not depend on the quench duration St. This result led us to conjecture that this 
late time behaviour agrees with that in an instantaneous quench. 
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In this section, we return to this issue, first by reviewing the d = 3 result and 
then by considering late time behaviour for higher dimensions. We will show that 
the agreement between smooth and instantaneous quenches found in d = 3 does not 
generally occur in higher dimensions. 


5.1 Review of d = 3 

The starting point will be to consider the correlator in eq. (2.17) for instantaneous 
quench, and evaluate this expression at coincident points in space and time, i.e., x = x', 
t = t'. Of course, this expectation value is divergent in the UV, so it must be regulated. 
In [2] we showed how to carry out the regularization in detail, but for now it will be 
enough to compute the difference between the quenched expectation value and that for 
a fixed mass m to produce a finite result. After subtracting, we get the finite difference 


m 


(0 '} quench (0 } fixed . 


dk 


dvr J k\/k‘^ + m? 


sin^ kt 


(5.1) 


for d = 3. Interestingly, this expression can be integrated analytically and the solution 
expressed in terms of generalized hypergeometric functions. 


(0 ^quench (0 ) fixed 


rri^t^ 


dvr \2t 


IT 


— 1F2 ( 1 , 


2-^3 


1 , 1 ; -, -, 2 ; 
2 2 


(5.2) 

The complete analysis of this solution can be found in [2], but let us just say here 
that the expectation value begins by growing linearly when mt ^ 1 (but still t/St 3> 1) 
but then for very late times, i.e., mt 3> 1, the expectation value keeps increasing 
but now only at a logarithmic rate, i.e., {(f>'^)ren ~ log(mf). As we show in [2], the 
instantaneous quench and the smooth quench calculation coincide for d = 3, basically 
because the integrand for the smooth quench decays fast enough, in a way that the 
approximation of u6t -C 1 continues to be valid. Recall from the discussion at the end 
of section 2, that it is possible to obtain the instantaneous quench expectation value 
starting from the smooth quench and taking both the late time limit, i.e., t/6t S> 1, 
and the low energy limit, i.e., u6t -C 1, for every u in the problem. However, we 
generally need to integrate momentum k up to inhnity with hxed 5t, so usually this 
approximation will break down for large enough k (remember that ujout = k m the 
quench to the critical point). In the special case of d = 3, though, the integrand decays 
in a way that only the low momentum modes contribute and then the approximation 
is reasonable. In the next subsection, we will show that this is a particular effect of the 
three-dimensional case and does not hold in higher dimensional spacetimes. 
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5.2 Higher dimensions 


The main problem that arises in taking the late time limit in higher dimensions is 
that the expectation value for the instantaneous quench cannot be regulated by simply 
subtracting the hxed mass expectation value for d > 3. Moreover we will show that it 
cannot be regulated using the usual counterterms found in [1, 2], This fact will be taken 
as a hint to argue that in fact, the low energy approximation is not valid in evaluating 
the late time expectation value of 0^ in higher dimensions. Then, in order to get the 
expectation value for the smooth quench, what one should really do is to fully evaluate 
eq. (2.13) in the limit where t/St 3> 1. Of course, without taking the extra low energy 
approximation, it will be impossible to recover the instantaneous quench answer, that 
it will turn out to be UV divergent for d > 7 and then, inhnitely different from the UV 
hnite result for the smooth quench. 

So, let us start by evaluating the bare expectation value for 0^ in the case of an 
instantaneous quench. In this case we have, 

i^'^) = f = (P + sin^(H)) . (5.3) 

^ 2(27r)'^-i J ^ ’ 


To explore the UV behaviour, we expand the expression above for large k and for up 
to d = 9, the results can be summarized as 


( 0 ^) = — / dkk 
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Of course, the terms appearing in the hrst line are those same divergent terms that 
we expect from the constant mass case. 


(0 ) fixed 




-2 


^/w+ 


zdk. 


m^ 


(5.5) 


But in eq. (5.4), we also have divergent terms in the second line proportional to 
sin^(A;t) that do not correspond to any physical counterterm contributions found in 
eq. (3.11). In fact, as we are interested in the long time behaviour of the expectation 
value and given that the mass prohle (and its time derivatives) decay exponentially 
in time, the only remaining physical counterterm in this limit should be the mass 
independent term of eq. (3.11), i.e., k^~^. 

In evaluating {(l)‘^)ren, we integrate over all momenta, but the integral is divergent 
(even after taking into account the physical counterterms). The reason for this be¬ 
haviour is that in higher dimensions the approximations of eq. (2.18) are not longer 
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valid. One way to realize this is to compare it with the expression before that approx¬ 
imation. Recall that this is given by eq. (2.13), which after some manipulation, in the 
late time limit becomes. 



where and (3^ are given by eq. (2.10). Now we wish to compare this integral with 
the result for an instantaneous quench 


= IT.' ^ J(Ik sin“(*:()) -k‘‘ . 

(5.7) 

For this comparison, we start by examining the integrands <F^(fc). For d = 5, this is 
done in £g. 8, where we choose 6t = 10“^ and we evaluate the expression at f = 10. 
However, our results generalize to the full range of values where the approximation 
of late times is valid. What we see is interesting: if we focus on small momenta, as 
shown in £g. 8(a), both the approximate and the full integrands agree. They both show 
a highly oscillating behaviour, that seems to continue to larger momenta and which 
would make both integrals diverge if it did so. However, what we see in £g. 8(b) is 
that in fact this behaviour does not continue for very large k in the case of eq. (5.6). 
It can be seen that <h^(/c) decays to zero for large momentum in the case of the full 
integral, as required to produce a UV hnite result. Instead, the approximate integrand 
continues to oscillate and so produces a UV divergent integral. In fig. (8), we see that 
the two integrands differ substantially for k > 1000, i.e., for k5t > 1. Of course, the 
approximation of Uout^t -C 1 is no longer valid in this regime and hence it is natural 
to expect that they should differ there. The approximations of eq. (2.18) are not 
valid to obtain the correct late time limit of (0^) in higher dimensions and hence the 
expectation value does not match that after an instantaneous quench. To complete 
this analysis we show that the same happens for d = 7 in £g. 9, where the oscillatory 
behaviour is even increased by a power law divergence in the approximate integrand. 
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(a) for small momenta. 


(b) for large momenta. 


Figure 8. (Colour online) Analysis of the approximation of low energies and late times in 
d = 5. We show results for dt = 10“^, t = 10 and m = 1. In the first figure, we show that 
in both eqs. (5.6) and (5.7) coincide when evaluated for small k. However in the second 
plot, we see that the full integrand (blue) decays for large k while the approximate solution 
(red) keeps oscillating. The second plot looks fully painted because of the highly oscillatory 
nature of the functions. 



(a) $^(A:) for small momenta. 


(b) $^(fc) for large momenta. 


Figure 9. (Colour online) Analysis of the approximation of low energies and late times in 
d = 7. We show results for dt = 10“^, t = 10 and m = 1. In the first figure, we show that 
<h^(fe) of both eqs. (5.6) and (5.7) coincide when evaluated for small k. However in the second 
plot, we see that the full integrand (blue) decays for large k while the approximate solution 
(red) keeps oscillating. The second plot looks fully painted because of the highly oscillatory 
nature of the functions. 


To conclude this section, let us summarize our findings with regards to the late 
time limit after the quench. We have considered two different protocols to quench a 
scalar held: the instantaneous one, where we suddenly start evolving an eigenstate 
of the massive case with the massless Hamiltonian; and the smooth one, where we 
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continuously evolve the mass of the scalar held from the massive case to the massless 
in a time scale of 6t. Now one would think that these different protocols must give 
different answers in the early time evolution,® but that in the limit of 6t going to zero 
and for late times, we should obtain similar results. In [2] we showed that effectively, 
if we take the limit of t/6t ^ 1 and also u6t -C 1 for every u, we can reproduce the 
instantaneous quench result from that for the smooth quench and so in principle, we 
might expect the same late time behaviour in both cases. This allowed us to identify 
the interesting logarithmic growth behaviour of the scalar held for late times in d = 3. 
In a self consistent way, we showed that both of these approximations were reasonable 
in d = 3 and so, the late time behaviour for the smooth and the instantaneous quench 
agreed. However, we found that this agreement in d = 3 was fortuitous because when we 
tried to repeat the analysis in higher dimensions, we found that the approximations of 
eq. (2.18) are no longer valid. That is, higher momenta (and hence, higher frequencies) 
contribute signihcantly to the expectation value of in higher dimensions and cannot 
be neglected. Hence for d > 4, the instantaneous quench gives a result for (0^) at 
late times which is inhnitely different from the smooth quench result. In particular, 
the smooth quench gives a hnite late time limit for (0^) as df —)■ 0 [2], while the 
corresponding result appears to diverge in an instantaneous quench. 


5.3 Regulated instantaneous quench 

It is interesting to note that the integrand in eq. (5.7) for the instantaneous quench 
in higher dimensions does not decay to zero for large momentum. Instead, it seems to 
show a rapid oscillatory behaviour around zero, as shown in hgs. (8) and (9). So one 
may think that even though the amplitude is diverging, the positive and negative part 
are cancelling in every period and so, in some sense, we may be able to recover a hnite 
result from these integrals. In fact, we are inspired here by the way that the hxed mass 
correlators were regulated in Appendix A. 

Hence let us go back to our instantaneous quench results. For d = 5, the behaviour 
of the expectation value of for large k can be extracted from eq. (5.4). This gives, 

(0^) -)■ — [ dk ({kt)\ . (5.8) 

fe->-oo as J \ 2 / 


The hrst divergence, proportional to k'^ ® is our usual counterterm that we will subtract. 
But note, then, that the term proportional to sin^ can be re-written to yield 
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®In our case, we obtain this interesting universal scaling near t = 0. 
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and this is one of the integrals that we now know how to regulate (just put a = 0 and 
X = 2tm eq. (A.4). Of course, this is just showing the large k behaviour of the integral. 
In order to get the full answer we should include all momenta and this, unfortunately, 
can only be done numerically. But, in principle, since we know how the integral behaves 
for large k, we should be able to get a hnite result for our integral as well by using 
this new regulator. The results are shown in £g. 10, where we evaluate numerically the 
instantaneous quench solution of eq. (5.7) for d = 5 and adding a regulator exp(—afc), 

{4)“^{a))instant = (^7^ j dk + m^sin^(fct)) - k^^ exp{-ak). (5.10) 



Figure 10. (Colour online) Expectation value of 47 in the instantaneous quench as a function 
of the regulator parameter a for d = 5 and mt = 10. The values of a go from 10“^ to 10“^ 
in a logarithmic scale. The solid line represents the value for the smooth quench at mt = 10 
and 5t = 1/20. 


We evaluate the expectation value at late times, mt = 10, and compare it with the 
smooth quench also at late times, where we are using St = 1/20. We hnd that as we 
take a to zero, the regulated integral for the instantaneous quench approaches some 
hnite value, showing that the integral converges. However, this value differs from that 
for the smooth quench in a relative amount by 


(0 y 0))instant (0 )smooth 


0.1670 - 0.1490 

smooth 


0.1490 


~ 0 . 121 . 


(5.11) 


So even if we found a way to make sense of the divergent integral in d = 5, the result 
does not quite coincide with the smooth quench result, where no approximation is made 
(apart from the late time limit approximation). 


28 













Note that the relative difference in eq. (5.11) is of order m5t. To quantify this 
difference more precisely we compute the expectation value in the instantaneous quench 
case with a = 10“^ and mt = 10. We then vary 6t in the case of the smooth quench and 
compute the relative error between the two results. The outcome is shown in £g. 11. 
We see that for ~ 1, the two protocols give very different answers. But this is 
expected because a large 6t means going into the adiabatic regime and this need not to 
agree with the rapid quench even at late times. Instead, as we decrease St, we see that 
the relative difference between the two approaches also diminishes and in fact, when St 
is of order 10 “^, the relative error is also of that order of magnitude. 

Naively, one may think that it is possible to understand this behaviour by expand¬ 
ing the expectation value in powers of St. To do that we start with the smooth quench 
integral in eq. (5.6). We then expand the Bogoliubov coefficients ag and j3^ for small 
St and compute the integrand to lowest orders in St. This results in 

1 f 

smooth = — dk {{k^ + rn^ sin(kt)^) + (5.12) 

0's J yk^ -\- 

+St^ ^ 12 ~ sin^(H)) -I- 0{St^) 

Let us analyse this last expression. The first term is independent of St and one can easily 
see that it exactly matches the expression for the instantaneous quench in eq. (5.7). 
Of course, this match was known implicitly, since expanding for -C 1 in the previous 
expression is actually expanding for uSt -C 1 and the agreement found above is the 
claim that we could reproduce the instantaneous result by taking the small frequency 
limit of the smooth quench. The next term in the St expansion appears at order St^. 
So, again naively, one might expect that smooth = {4^'^)instant + 'ySt'^ 0{St^), for 
some number 7 . However, if we look carefully at the integral that gives this correction 
at order St"^, we will see that it is in fact divergent for d > 5. We can try to regulate 
it by adding a regulator as in Appendix A, but in fact we will see that apart from 
the oscillating term (which can be regulated) there is an extra constant term in the 
integrand, i.e., — that will make the integral divergent as we integrate over 

k from 0 to 00 . 

This is closely related to the fact that if we try to £t the relative error by some 
power law expression in the region of small St — see red dashed line in £g. 11 — we 
hnd that the error does not scale as St^ but the exponent is rather close to 1.40. This 
is another sign that this naive expansion is somehow ill-dehned for d = 5. Again, we 
should say that behind this expansion it is assumed not only that St is small but that 
uSt <C 1, for every uj in the problem, and we already showed that the assumption is 
not valid for d > 5. 
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Figure 11. (Colour online) Relative error in the computation of {(fP') at late times using 
the smooth and the instantaneous quench as a funtion of 5t. Here we are regulating the 
instantaneous result with a = 10“^ and evaluating at time mt = 10. The red dashed line is 
showing the fit of the first points (from 6t = 1/1000 to 6t = 1/800) by a power law function 
of the type f{6t) = where a = 20.03 and a = 1.404. The inset zooms in the region 

where the fit was made showing perfect agreement between the points and the fit. 


The situation is even worse in higher dimensions, where the hrst term in the series, 
i.e., the order St^ term in eq. (5.12), fails to converge. As an example, we show what 
happens in d = 7. Even though in Eg. 9 the integrand appears to oscillate around zero, 
this is not the case. For large /c, we have, after subtracting the usual counterterms, 

(0^)ren-1 — f dk + (rfpk'^ - sin^{kt)\ 

k^oc (Jg J \ I 0 \ 2 J J 

~ — y dk — 2rrPk‘^) cos(2/ct) + . (5.13) 

So, even though the hrst two terms could be regulated using the above prescription, 
there is an extra constant term that cannot. Adding the regulator to the term propor¬ 
tional to m^/8, we will get a result which diverges as a —>■ 0, 

4 poo 4 

tI (5.14) 

and so, the limit of a —)■ 0 in this case will be nonsense. The same also happens in any 
higher number of dimensions. In fact, the only reason why this worked in d = 5 was 
because the divergent terms with the sin^ term, exactly matched the terms without 
the sin^ term, in a way that made the whole integrand to oscillate around zero. But 
this, as shown above, does not happen in general and so, the instantaneous quench 
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approximation in higher dimensions gives a value for ( 0 ^) which differs by an inhnite 
amount from the smooth quench result, even in the late time limit. 

6 The energy density at late times 

In [1, 2], it was argued that the renormalized energy density for a quench satisfying eq. 
(1.2) obeys a scaling relation (1.3), 

. ( 6 . 1 ) 

This result is consistent with the scaling of the expectation value of the operator (1.4) 
since it satishes the Ward identity 

^ = -dt\{t){0^)ren. (6.2) 

In the corresponding scaling relation (1.4) for the quenched operator, {0^)ren is uiea- 
sured earlier than or soon after the end of the quench. However, the energy scaling (6.1) 
will be valid for arbitrarily late times since the energy is injected into the system only 
during the quench. The equation (6.1) gives the 6t dependence. The energy density 
itself will have additional hnite pieces, which would be subdominant for A > d/2, but 
in fact give the dominant contribution for A < d/2. 

In this section, we will concentrate on the energy density for the free bosonic 
held with the mass prohle (2.2) at asymptotically late times and calculate a UV hnite 
quantity: the excess energy above the ground state energy of the system with the value 
of the coupling at asymptotically late times. We will perform the —)■ 0 limit and 
compare the results with that of an instantaneous quench. 

In terms of the “in” modes the energy density is given by 

^ = \f (^Try-i (6.3) 

Since we are interested in the behavior of this quantity at late times, it is convenient 
to express this in terms of the “out” modes. In terms of the Bogoliubov coefficients 
f^k dehned in (2.9) and (2.10) this becomes 

^ = \ j ^2 -kY-^ + + + 

+ + m^{t) - uU{a^l3lv^v_j: + apj^vp/^)). (6.4) 
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In deriving this expression, we have used the fact that and depend only on \k\. 
At t = oo, the second line of eq. (6.4) vanishes. Using the relation \cx^\^ — = 1 

and the asymptotic form of the out modes at f — )■ cxd, i.e., 




y/^A^Jowt 


^i(k-x—uJoutt) 


(6.5) 


one gets 

■ ( 6 . 6 ) 

However, the ground state energy of the system with the hnal value of the mass is given 
by 


^ _l f d‘^-^k 

Oground ~ 2 J ' 

Therefore the excess energy over the hnal ground state is given by 


(6.7) 


AS 


d'^-^k 

{271^-^ 


1/3^1 


( 6 . 8 ) 


Using the explicit form of the Bogoliubov coefficients in eq. (2.10) and integrating over 
the angles, we arrive at the hnal expression 


AS = 


Sld- 


d-2 


( 2 ir: 


d-1 


dk k 


d-2 


^out 


sinh^(7ra;(5t) 


sinh(7ra;j„(5t) sinh(7ra;o„t(5f) 


(6.9) 


For the mass prohle (2.2), the integral in eq. (6.9) is hnite for any hnite dt. In fact, 
for small fc, the integrand approaches while for large k, it becomes 

(^d-3)^4(^^)2g-27rfc(5t_ Heuce the integral above is convergent both in the IR and UV 
for any physical d >2. 

To analyze the small A = m6t limit of eq. (6.9), let us hrst scale out a power of 6t 
to write AS = {St)-‘^Ii{X) where 


7i(A) 



sinh^[|(^g2 y2 _ 

sinh(7riyg2 _l_ sinh(7rg) 


( 6 . 10 ) 


Clearly /i(0) = 0. However the small A dependence is diherent for diherent dimensions. 
It turns out that 


/i(A) ~ A'^ c/ = 2,3, 

~ A^ d>4. (6.11) 

The above behavior was determined from a direct numerical evaluation of the integral 
in eq. (6.10) and htting the results shown in hg. 12. 
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Figure 12. (Colour online) The integral /i(A) in eq. (6.10) as a function of A on a logarithmic 
scale. The data are for d = 2,3,4, 5,6, 7 from top to bottom. The solid lines correspond to 
/i(A) ~ A'^ for d = 2,3 and /i(A) ~ A^ for d = 4,5, 6, 7. Note that the fit in d = 4 is not as 
good as the others, which suggests that there are log corrections to the leading behavior for 
d = 4. 

This means that for d = 2, 3, the excess energy has a smooth finite small 6t limit 
with A£ ~ m'^. The leading answer is exactly the same as the energy excess for the 
instantaneous quench, which can be read off easily from the corresponding Bogoliubov 
coefficients (2.16) (for the mass profile ( 2 . 2 )) 



( 6 . 12 ) 


This quantity is finite for d = 2,3 and hence we have 



(6.13) 


To estimate the corrections to this leading small A behavior consider the difference of 
the excess energy to the excess energy due to an instantaneous quench, 


SS = A£- A£ 


•instant 


St-'^ (Ji(A) - (6.14) 


•instant 


A numerical evaluation of this quantity shows 


5 £ ~ , 


(6.15) 
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Figure 13. (Colour online) The absolute value of / 2 (A) = m~^ 5£ as a function of X = m6t on 
a logarithmic scale. The data are for d = 2,3 from bottom to top. The solid lines correspond 
to the curves / 2 (A) ~ A^ for d = 2 and A for d = 3. Those curves correspond to the best fit 
of the numerical data in the region 10“^ < A < 10“^. 


both for d = 2 and 3. This behaviour is shown in £g. 13, which is a log-log plot of the 
quantity / 2 (A). Clearly, we have / 2 (A) ~ A^ for d = 2 and / 2 (A) ~ A for d = 3, which 
leads to the scaling in eq. (6.15). 

On the other hand, for d = 4,5,..., one recovers the desired scaling from the 
leading behaviour, i.e., 

AS ^ XUt-‘^ = , (6.16) 

with logarithmic corrections for even dimensions. Thus for these dimensions, the energy 
density diverges in the m6t —)■ 0 limit. At the same time, the energy density for an 
instantaneous quench diverges in the UV for these dimensions. Let us emphasize this 
point once again: our results show that the energy density is UV hnite after a smooth 
quench; however, in the limit mSt —?• 0, that energy density diverges, just as in the 
instantaneous quench, where the divergence is in the UV. 

These results also indicate that for d > 4, the excess energy is in fact given by linear 
response. This is in accord with the computation of the renormalized expectation value 
of and the renormalized energy density for d > 4, as described in [1, 2]. Indeed if 
we extract the leading piece in a small A expansion by expanding the integrand in eq. 
(6.10) we get the expression 


/i(A) ~ 



Q 


d-3 


(sinh(7rg))2 


(6.17) 
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This expansion makes sense when the above integral is convergent. However the integral 
is IR divergent for d = 2,3,4. The latter explains is why /i(A) does not have an 
expansion in terms of in these dimensions. 


7 Excess Energy for General Theories 

The discnssion in section 6 snggests that the scaling of the excess energy shonld be 
a property of a general interacting held theory. In this section, we argne that this is 
indeed trne. 

Consider an action 



(7.1) 


where Soft is a conformal held theory action. The fnnction X{t) is of the form 



Alternatively, we may write X{t) = Aq + 6XF{t/6t) if we specify F{y < 0) = 0 and 
F{y > 1) = 1. We leave the details of the fnnction F{y) dnring the transition [i.e., 
0 < 1/ < 1) nnspecihed other than that the maximnm is hnite with F^ax > 1- Farther, 
this prohle may dip below zero by some hnite amonnt and so we specify the minimnm 
as Fjnin < 0. Implicitly, we are also assnming that the prohle is smooth. 

The system is prepared in the gronnd state of the initial action. Let us evaluate 
the total energy density £{t) at some time t > St in a perturbation expansion in SX. 
To quadratic order in SX, this expression is given by 


S{t) =£^-SX F{t/St) (0, Ao|0(0,0)|0, Aq) 



^0 



0 Jo 

, i/ 


dt' I dt" / d'^-^x’ F{t'/St)F{t"/St) {t),Xo\0{x',t')0{Q,t")\t),XQ) 



-£qSX^ I dt' I dt" d'^-^x'F{t'/St) F{t"/St) {0,Xo\0{0,t")0{x',t')\0,Xo) 


dt' / dt" / d‘^-^x'F{t'/St)F{t"/St){0,Xo\0{x',t')HoO{6,t")\0,Xo) 


+0{SX ^), 


(7.2) 


where Hq is the initial Hamiltonian, |0, Aq) is the ground state of the initial Hamilto¬ 
nian, and So denotes the initial ground state energy density. Here we have used space 
translation invariance as well as the fact that one point functions in the initial ground 
state are constants in both space and time. On the other hand, the expectation value 
of the operator O is, to order 6X, given by 

(0,A„|C'(f,i)|0,A„) = {0.A„|O(0,0)|0.A„)-! j d‘‘-'x j dt'Mt') (0. A„||C'(f, i), 0(0, i')l|0, A„>. 

(7.3) 

Using this, it is straightforward to verify that the Ward identity 

^ = -^{0,A„|O(f,()|0,A„), (7.4) 

is satished.® 

We are interested in evaluating eq. (7.2) at late times. However since the coupling 
is a constant for t > 6t, the energy density at inhnitely late times is exactly the same as 
the energy density at t = St. The ground state energy density of the hnal Hamiltonian 
is given by the standard expression 


Sf — Sq — SA (0, Ao|0(0, 0)|0, Aq) 




E 


l(0,Aolfd^ 0)|n, Ao)p 

So — Sn 


(7.5) 


where Sn denote the energy densities of the excited states \n, Aq) of the initial Hamil¬ 
tonian and V is the volume of the system. It is clear from eqs. (7.2) and (7.5) that 
the excess energy density AS = S{St) — Sf starts at 0{SX^). Moreover we expect that 
the UV divergences in S{t) for t > St are cancelled by those in Sf. This expectation 
comes from the following fact: as seen in the previous sections, and in [1, 2], the UV 
divergent terms depend on both X{t) and its time derivatives. However for t > St, the 
coupling is constant and these time derivatives vanish. Therefore, the UV divergent 
terms should be those of the constant coupling interacting theory. While this is explicit 
in the free held theory considered in section 6, we do not have an explicit proof for 
general interacting theories but it stands as a reasonable expectation. If the resulting 
expression for AS is also IR hnite, it is determined to this order entirely by dimensional 
analysis, 

AS SXSt<^-^^ . (7.6) 

Further, as discussed in [1, 2], the corrections to this result would be a power series in 
the dimensionless coupling g = SXSt^~^, which is small in the fast quench limit. 

® The time derivatives of the second and third lines in eq. (7.2) cancel the time derivatives which 
act on the upper limit of integration in the third term of the first line. 
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This argument will fail if the integrals involved in AT are IR divergent. This can 
be seen to happen when 2A > d, as we have explicitly seen for the free held theory for 
d = 2 and 3. 

8 Discussion 

The aim of this paper is to establish a precise relation between the smooth fast quench 
and the instantaneous quench. These are the most common quench protocols discussed 
in high energy theory and condensed matter physics literature, respectively. Naive 
reasoning would say that if one considers the evolution at very late times (with respect 
to the quench rate 6t), then both protocols should give the same result since 6t would 
be negligible. However, our results in this paper suggest that they may or may not give 
the same answers depending on a variety of factors such as the spacetime dimension, 
the scaling dimension of the quenched operator and how much time after the quench 
is considered. In this paper, we computed spatial correlators, local expectation values 
and the energy density. In this section, we summarize our results, making precise 
statements on when the abrupt approximation makes sense. We also discuss different 
procedures to regulate the energy density. The bulk of our comments below relate to 
the explicit calculations performed in the free field theory. However, as we will discuss 
at the end of this section, we expect that these conclusions hold for generic interacting 
theories for smooth fast quenches as dehned in eq. (1.2). 

Spatial Correlators 

The study of spatial two-point correlation functions is interesting because they are UV- 
hnite quantities that introduce a new scale to the problem, i.e., the spatial separation 
r. The behaviour of this object is very different depending on the time at which we 
compute it. We summarize the results here for early times, i.e., t/6t ~ 0(1), and for 
very late times, i.e., t/6t ^ 1 and mt S> 1. 

• At early times, we can distinguish between three different regimes. 

When m6t > 1, independently of r, we are in the slow quench regime, so we cannot 
compare with our previous results. However, there should be some signatures of 
universal behaviour corresponding to the Kibble-Zurek scaling when quenching 
through a critical point. We leave this appealing point for further research in the 
future [16]. 

The most interesting feature appears when r < St < 1/m. This means that 
the quench is fast since mSt < 1 but the spatial separation dehnes the smallest 
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scale. Then we found that the same universal scaling that was reported in [1, 2] 
appears in the two-point spatial correlation function. This holds in any spacetime 
dimension and for arbitrary “early” times. We can think of this correlator as a 
version of (0^) regulated with point splitting and so, r plays the role of a short 
distance cut-off — see below. 

In general, we would continue to decrease 6t towards the instantaneous limit 
in which 6t —)■ 0. However, in this correlator, we are limited by the distance 
separation r. In fact, the correlator saturates as 6t gets of order r and then the 
result becomes independent of 6t. In all this analysis we were able to take the 
UV cut-off to inhnity but we expect a similar behaviour when working in theories 
with a hnite cut-off, with r~^ playing that role here. 

• At late times, the results depend on the spacetime dimensions and on the sep¬ 
aration distance. For long distances, in any dimensions, the correlator for the 
instantaneous quench and the smooth quench coincide. As the separation be¬ 
comes smaller, the behaviour is different depending on the spacetime dimensions: 
for d = 3, the smooth and the instantaneous correlator continue to coincide as 
r —)■ 0; for d = 5, there appears a small hnite difference that goes to zero as 6t 
goes to zero; hnally, for d = 7 the two correlators differ by an inhnite amount as 
r —)■ 0. We expect the latter behaviour extends in higher dimensions. 

Expectation value of (fp' 

We showed that the short distance expansion of the correlator is in one-to-one corre¬ 
spondence with the counterterms needed to regulate the bare expectation value of (fP. 
Then, it should not be a surprise that the behaviour of the {(fP) at late times is very 
similar to that of the spatial two-point correlator at late times but with small spa¬ 
tial separation. In fact, we showed that for d = 3, the smooth and the instantaneous 
quench give the same answer. When evaluated for d = 5, the smooth quench differs 
from the instantaneous quench but only in a hnite amount that is of order St. In higher 
dimensions, however, it is impossible to regulate the expectation value of (jP in the case 
of the instantaneous quench. The smooth quench, in contrast, has a smooth hnite limit 
as —?■ 0 and so, the two approaches yield inhnitely diherent results. 

Regulating the energy density 

Both in the present and previous [1, 2] works, we worked in a framework where UV- 
divergent quantities where regulated by adding suitable counterterms. Of course, we 
showed how to construct those counterterms and how they yield hnite values for quan¬ 
tities such as the expectation value of ipP) and the energy density. The way in which 
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such a subtraction is done is much in the spirit of how regularization works in AdS/CFT 
through what is known as holographic renormalization [17-19]. There we add an extra 
counterterm (boundary) action to the usual gravitational action to get hnite expec¬ 
tation values. This is how, for instance, expectation values for holographic quantum 
quenches are regulated in [12], which served as a motivation to our studies. Our ap¬ 
proach is also reminiscent of the way held theories in curved spacetimes are regulated. 

However, a number of things about this procedure may appear strange to a typical 
held theorist. In particular, the fact that our counterterms hrst, depend on time and 
second, some terms depend on time derivatives of the quenched coupling. In this 
section, we would like to go back to this procedure and compare it with other candidates. 

To summarize, we dehne a renormalized energy density by subtracting counterterm 
contributions from a bare energy density (and taking the cut-oh to inhnity). Basically, 
(i^)ren = ^quench “ ^ct where Squench and Set Separately diverge as the cut-oh goes to 
inhnity but {S)ren is hnite. We also showed in [2] that with this dehnition of renormal¬ 
ized energy density (and an equivalent for the scalar held), we satisfy the Ward identity 
in any spacetime dimension and for any quench protocol. 

A second, perhaps more standard approach, would be to recognize the divergences 
as coming from the zero-point energy for the scalar held. Each momentum mode 
behaves as a single harmonic oscillator and then if we sum all the zero-point energies, 
i.e., we get a divergent quantity. This is what we call Sfi^ed, that is, the energy 

density for a scalar held of hxed mass at any instant of time. Again, one would naively 
say that we should get a UV hnite value if we subtract Squench — Sfixed- One nice thing 
about this is that if we go to very early or very late times where the mass is constant, 
there is a ground state and it has precisely zero energy density and any other state has 
a positive energy density. Even if this procedure works for low dimensional spacetimes, 
we showed that is not enough in higher dimension. In fact, for d = 3 and 5, we have 

Sfixed Set 0 ^ ^ ) (® 0 ) 

while in d = 7, 

s fixed - Set (Xm^ + 171^ A, (8.2) 

where A is some energy UV cut-oh. So hrst thing to note is that for d < 7 (actually, 
for d < 6) both the counterterm energy density and the hxed energy density only 
diher by a hnite amount, so if one is sufficient to regulate the theory then, so is the 
other. Moreover, as most of our study corresponds to the fast quench regime where 
mSt 1, this hnite amount would be negligible compared to the scaling with St and 
so, the conclusion will be unchanged using either approach. ' Let us also note that the 

^See, however, the discussion on the reverse quench in d = 3. 
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subtraction of S fixed also satisfies the Ward identity. This is easy to see as the only 
time dependence is on the mass, so 

diEfiad = dt j dkk'‘~'‘fk^ + j dk -^=T== j d,nd(t), 

(8.3) 

but the term in parentheses in the hnal expression is just the bare expectation value 
of 0^ with an instantaneous mass m{t) — compare to eq. (5.5)— , which gives exactly 
the Ward identity in eq. (2.12). 

The situation is completely different in higher dimensions, though. The energy 
density of the quench has more divergences than those appearing in £ fixed- These are 
proportional to time derivatives of the quenched coupling. If we suppose that all the 
quench happens within a time of order 6t of, say, t = 0, then these terms will not 
affect what happens at very early and very late times, so it is possible to compute 
£quench “ £ fixed In those regimes. However, if we want to follow the evolution through 
the actual quench, £ quench— £ fixed is just divergent and we do not have a hnite observable 
in the middle of the process. This is the main reason why £xen is a better measure of 
what is going on during the quench, i.e., because it allows us to compute the energy 
density at any time in any spacetime dimensions. 

All of these situations are depicted in £g. 14, where to show clearly what is going 
on we set as the zero of energy density with £et in the hrst row and with £ fixed in the 
second. For d = 3 and 5, both approaches are valid and we see that the only difference 
is on some small quantity proportional to m. We showed both d = 3 and d = 5 because 
£fixed and £ct differ in each case by a different amount. While in d = 3, the difference 
is negative, in d = 5, it is positive. In d = 7, however, we can follow £ren but not the 
other one. Then, it is clear now that the naive intuition is wrong or at least is not 
complete. 
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Figure 14. (Colour online) Schematic description of the evolution of the energy density after 
a quench as a function of time depending on the regularization scheme. In the first row we 
exemplify the counterterm subtraction that allows us to follow the evolution for any time, 
even during the quench itself, characterized by the time scale 6t. Taking the counterterm 
energy density as the zero of energy density, at very early times, the starting energy density 
is negative for d = Ak + 3 and positive for d = Ak + 1, with k > 0 integer. In the second row 
we take the hxed energy density as the zero of energy density. In this case, all the quench 
energies are greater than this ground zero energy density. The reason is that the fixed mass 
energy density is the energy density for the scalar with a fixed mass and no quench, so any 
energy inserted during the quench will give a positive additional contribution. Note that 
in low dimensional spacetimes, the fixed energy density differs from the counterterm energy 
density but always by a finite amount proportional to a power of the mass at that instant of 
time, while the quench energy density usually scales also with 5t, so our universal scalings, 
i.e., 6£ ~ 6}? will appear in any case. The same will happen in greater dimensions 

at very early and very late times. However, for d > 6 during the quench, it is not sufficient to 
subtract the fixed energy density, as there are extra UV divergences in the quenched energy 
density that are proportional to time derivatives of the mass. In this case, subtracting the 
fixed energy density is sufficient to compute the energy density at very early or very late times 
but not during the middle of the quench. This is depicted in the plot in the bottom-right 
corner for d = 7 but is a general feature of higher dimensional spacetimes. In the same way, 
in the upper-right plot we cannot sketch £fixed — £ct during the quench as they differ by an 
infinite amount. 
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The next step would be to think whether there is some other way to regulate the 
energy density in higher dimensions. A possible answer already appeared in [1, 2] while 
not completely emphasized. We know that the counterterms come from an adiabatic 
expansion. At zeroth order, this adiabatic expansion gives just the hxed energy density 
that corresponds to doing the quench inhnitely slowly so that, at each instant of time, 
the energy density is just the energy density needed for the scalar held to have that 
particular mass m(t). What we do to get the counterterms in just to expand the 
adiabatic expansion for large momentum and then extract the divergent pieces. As 
mentioned, at zeroth order, the hxed energy density dihers from the counterterm energy 
density by a hnite piece proportional to But we also know that for higher 

dimensional spacetimes we need to go to higher orders in the adiabatic expansion to 
capture the divergences involving time derivatives of the mass. So an idea to generalize 
the hxed energy density subtraction to higher dimensions would be that, instead of 
subtracting just the counterterms, to subtract the full energy density in the adiabatic 
expansion to that order. This would correspond to the energy of a slow quench but 
going beyond the zeroth order. 

To be more explicit, in the adiabatic expansion shown in [2] we dehned incoming 
modes of the form 

Wfc = ^7=^exp (ifc-f-i/ VLk{t')dt^ . (8.4) 


Then we found that 




1 

Auk 




+ 0 




(8.5) 


where So the zeroth order term in is the hxed mode energy density 

but this is then corrected with a second term that is second order in time derivatives 
of the mass and so forth. 

For d = 7, we showed that it is enough to expand the energy density to that order. 
So instead of regulating the energy density with the counterterm energy density what 
we can do is to subtract, at any t/5t, the whole energy density coming from that second 
order expansion. This will include, of course, the necessary terms to cancel all the UV 
divergences but it will probably introduce some extra hnite terms, in analogy to the 
extra hnite piece that the hxed energy density has with respect to the counterterm 
energy density in lower dimensions. Let us add that this subtraction is also consistent 
with the Ward identity.^ 

®To see this one should use the Uj: modes of eq. (8.4) to compute the energy density for the scalar 
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In all, we have two different consistent ways of regulating the energy density in 
our scalar held quenches. It would be interesting to consider whether there is some 
analogy of these two methods in holography. Our counterterm subtraction is clearly 


the counterpart of the holographic renormalization approach. But what would be 


the equivalent of subtracting the hxed mass energy density in an holographic setup? 
Well, it seems reminiscent of the old method of background subtraction in the early 
days of the AdS/CFT correspondence (see, for instance, [20]). Usually, divergences 
in holography appear as we take limits toward the boundary due to the divergent 
nature of pure anti-de Sitter spacetime as we take the radial coordinate towards the 
boundary. So, the hrst idea in holography, which was inherited from early semi-classical 
calculations in quantum gravity, to get a hnite renormalized quantity for some excited 
state was to subtract that same quantity but in the vacuum state, i.e., in pure AdS. 
Then both quantities will be divergent but their subtraction would be hnite and this is 
quite analogous to our hxed energy density subtraction. Later on, this procedure was 
replaced by the more rigorous method of holographic renormalization. 

There is an interesting ehect in the use of these two diherent approaches in d = 3. 
In this case, the scaling is special because instead of giving a diverging behaviour as 
5t —)■ 0, it gives a vanishing one. We analysed this case in [2], concluding that actually 
the energy density produced was given by 6{£)ren = where m is the initial mass. 
The interesting thing was that then, if we do the reverse quench where the initial mass 
is zero, it appears that the work done by the quench is zero! However, this was an 
artifact of using the counterterm energy density to regulate the expectation value. If 
we use the hxed energy density then we will hnd that the energy density starts from 
zero at early times and it goes to some hnite value, giving some non-zero hnite work in 
the process, as depicted in hg. 15. 


field, 



( 8 . 6 ) 


Upon taking the time derivative of that expression, one should get two different terms. First, a term 
proportional to 



(8.7) 


that perfectly vanishes since that is the equation which Qk should satisfy in order for the modes to 
satisfy the equations of motion (see [2]). The second term, however, is not vanishing and gives exactly 
the Ward identity. 
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Squench Set — Set Set — 0 


Sfixed Set 


(a) The energy density using the counterterm subtraction 



(b) The energy density using the fixed mass subtraction 


Figure 15. (Colour online) Schematic plots of the evolution of the energy density in three 
spacetime dimensions when 6t —)• 0. From the counterterm subtraction point of view, the 
work done is zero but from the fixed mass subtraction perspective, there is finite work done. 


Finally, we wish to consider one last method of obtaining a hnite energy density. 
In section 6, we considered the difference between two physical energies, i.e., Squench — 
Sground- Note that Sground, the ground-state energy density, should be something that 
we can easily dehne at very late and very early times. In particular, this looks like 
Squench — Sfixed at these early and late times. However, Sground is a real physical energy 
density of a particular state and it can be dehned in any renormalization scheme. So 
in this case, we do not need to make any reference to our choice of scheme because 
all of the divergences cancel in the difference of two physical energies. Of course, the 
drawback is that this can only be computed at very early or late times. 

Even dimensions 

Most of the explicit calculations presented in this paper refer to odd spacetime dimen¬ 
sions. However, most of the conclusions also hold for even dimensions. As pointed out 
in [1, 2], the only differences between even and odd dimensions are that there are addi¬ 
tional logarithmic UV divergences which must be regulated in even dimensions, and as 
a result, the renormalized expectation values have an extra logarithmic enhancement 
in the 5t scaling. For example, in even dimensions, the expectation value for 0^ in the 
fast smooth quench scales as 

(0^)ren ~ log/i5t, (8.8) 

where /i is a new renormalization scale introduced by the logarithmic counterterms. 
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The appearance of logarithmic counterterms adds an additional technical difficulty 
to the calculations presented in this paper but we do not expect that they would change 
the main results. In particular, the renormalized expectation values of have a smooth 
limit at late times as —)■ 0 in higher even or odd dimensions, while the analogous 
quantity diverges after an instantaneous quench. With regards to the energy density, 
we presented results for lower even dimensions in section 6, where we showed that the 
excess energy has a smooth limit in d = 2 as —?■ 0, which matches the instantaneous 
answer. In higher dimensions, the excess energy diverges as expected from the scaling 
of (0^)ren and the Ward identity (2.12). 

Lessons for interacting theories 

We end this discussion with a comment on the lessons of our work for general interacting 
theories. In [1, 2] we showed that the scaling form of renormalized quantities holds for 
general quantum held theories for fast quenches as dehned in eq. (1.1). It is natural 
to expect that the scaling for correlation functions found for the free theory in section 
4 would have an analogue in interacting theories as well. When the length scale in the 
correlator is small compared to the quench time, this correlator can be viewed as a 
point-split version of the operator which is used for the quench and in the fast quench 
limit, the arguments of [1, 2] then show that this quantity would scale in the expected 
fashion. 

In this paper, we found that the relationship between the fast limit of a smooth 
quench and an instantaneous quench is non trivial for free held theories in high dimen¬ 
sions. This again should generalize to interacting theories. What really led to the non 
trivial relation is the fact that in higher spacetime dimensions, the conformal dimen¬ 
sion of the quenched operator becomes large. Indeed, the scaling of the renormalized 
quenched operator O for general interacting theory together with the Ward identity 
shows that the renormalized energy density at late times behaves as and there¬ 

fore diverges as St ^ 0 for any d whenever A > d/2. This can be consistent with the 
results of an instantaneous quench only if the latter is UV divergent in this case. This 
fact should have non trivial consequences for the ability to express the state after a 
quench in terms of a boundary state as in [7] even in low spacetime dimensions when 
the conformal dimension of the quenched operator is large enough. 
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A Review of constant mass correlators 


In this appendix, we review the computation and behaviour of the (spatial) correlator 
for a massive free scalar held with a constant mass. For simplicity, we will focus on 
odd spacetime dimensions. Hence we are interested in computing the following spatial 
correlator, 

C'(f) = (0(f)0(O)) = —Axt f / ^ (A.l) 

^ ^ 2(27r)'^-i J + m2 ’ ^ ^ 

where m is simply a hxed constant (for all time). First, we can choose, without loss 
of generality, to place f along one particular axis using the rotational symmetry of the 
problem. Integrating out the transverse angular directions, then yields 


C(f) 


2(27r)'^-i 


k^-^dk 


dOsin 9e 


d—3 nAkr cos 0 


(A.2) 


where k = \k\ and r 
we hnd 


|rj. The integral over 9 can be done analytically and for odd d. 


C(r) 


1 

d-3 

a^r 2 


k^dk 


J d- 3 {kr) , 


(A.3) 


where Uc = 2 2 tt 2 and J d -3 is the Bessel function of order 

2 ^ ^ 

Now, to get the full answer for C'(r) we need to integrate over all k, so we note 
that for large k, the Bessel function behaves as times some linear combination 

of trigonometric functions (of kr). Hence the naive counting of the powers of k would 

d—4 

yield an overall factor of k^^ in the integrand above and hence one might conclude that 
the integral would diverge for any d {> 2). However, this factor provides a envelope 
for a rapidly oscillating function which tends to produce an added cancellation in the 
integral. Integrals of this form can be dehned with the following regulator: Insert an 
additional factor of the form exp (—a A;) to the desired integrand. The resulting (hnite) 
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answer now remains finite in the limit of a —)■ 0. In fact, one can show that this method 
works correctly for integrands that are a product of some power of k times an oscillatory 
function around zero. For example, one can show that. 


dk k'^ sm{xk + 6) = lim / dk k°‘sm{xk + 6) exp{—ak) 


a^O 


= sin + 1) + r(a + 1) x , (A.4) 


for any non-negative values of a and x. Applying this apprach to eq. (A.3) yields the 
following analytic answer 

( 7 ) ' K,_i{mr) , (A.5) 

where now is the Bessel K function. Of course, one can readily verify that this 
correlator satishes the Klein-Gordon equation as desired.® 

Given this expression (A.5), it is straightforward to examine various asymptotics 
of the correlator. In particular, considering the limit mr —)■ 0, we obtain 


acC{f) 


r(i-l)2" q_ mV _mV_ 

7Fr'^-2 V 2(d-4) 8(d-4)(d-6) ^ ^ 

r (1 - 

H- + 0{m'^r‘^). 


+ 


(A.6) 


Hence the correlator diverges as expected as r —)■ 0, i.e., the leading divergence goes as 
1/r'^-^. The above expansion reveals that the sub leading terms are all proportional to 
the mass in this limit. To have some concrete examples, we show: 


d = 3, 
d = 5, 
d = 7, 


(TcG(f) 

acC(f) 

acC{f) 


- m + Oirn^r ), 

r 

1 m? nfi . 

——+ 0 m r , 
r^ 2r 3 

;:a - TT-j + ^ - TIT + 0(mV). 


2r^ 8r 


15 


(A.7) 


We also note that the leading term in eq. (A.6) is, in fact, the exact answer for the 
massless correlator, i.e., 


C(fi = 


r(¥) 


^ i rd-2 

4:712 r 


for m = 0 . 


(A.8) 


^The correlator (A.5) applies for general spacelike separations if we replace the distance r by 


^j.2 _ ^ 2 ^ 
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The other interesting limit to consider is mr —)• oo. In this case, the Bessel function 
decays exponentially and we hnd 

d-3 

p rppf 2 

m = -^ (l + 0(l/(mr))) . (A.9) 

(Tc r 2 

We might note that the power of r in the leading term in eqs. (A.6) and (A.9) happens 
to coincide for d = 3 but otherwise they differ. 

To conclude this appendix, we emphasize the two main results: The hrst one is that 
naively the integrals above seem to be divergent, especially for high d. However, because 
the integrand is mainly oscillating around zero, they can be regulated as in eq. (A.4) to 
get a hnite result. The second lesson is that in the static case this correlator diverges 
as as shown in eqs. (A.6) and (A.7). We will take these facts into account when 

we analyse spatial correlators in the instantaneous and smooth quenches in section 3. 
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